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LINEAR AND NONLINEAR ANALYSIS OF 
ORBITAL TELESCOPE /SPACE SHUTTLE 
DYNAMICS AND CONTROL 

by 

S.M. Joshi^ 

SU>MARY 

During extended space shuttle orbital missions of the 1980 's, a number of 
solar-, stellar-, and Earth-viewing scientific experiments shall be performed. 
These experiments require a highly accurate instrument pointing system for 
providing fine pointing in the presence of crew-motion and other disturbances, 
and sensor noise. In view of this, an annular suspension and pointing (ASP) 
system, which makes use of a magnetically suspended vernier pointing assembly, 
is being developed for ultimate inplementation and flight testing, possibly during 
early space shuttle test flights. The objectives of this study are as follows: 

• Development of a detailed mathematical model of the space shuttle/ASP 
system. 

• Design of control laws in order to obtain the desired pointing 
performance. 

• Prediction of the statiscal pointing accuracies in the presence of 
stochastic disturbances such as crew-motion, and sensor and actuator noise. 

This report presents a part of the effort, consisting of linear analysis 
of the space shuttle/ ASP system which is currently in progress. The material 
presented includes derivation of a mathematical model and design of appropriate 
control laws. 


1 . iNTRODUCriGN 


During extended space shuttle orbital missions of the 1980*s a number of 
solar-, stellar-, and Earth-viewing scientific experiments shall be performed. 


^Research Assistant Professor in Electrical Engineering, Old Dominion University 
School of Engineering, Norfolk, Virginia 23508. 


These experiments will require a highly accurate instrument pointing system 
for providing fine pointing in the presence of crew-motion and other disturbances. 
With this objective in mind, a concept for such an apparatus, an annular 
suspension and pointing (ASP) system, was conceived (ref. 1). 

The annular suspension and pointing (ASP) system includes two assemblies 
with connecting interfaces, each assembly having a separate function. [See fig. 
1(a).] The first assembly is attached to the carrier vehicle and consists of 
an azimuth gimbal and an elevation gimbal which provide "coarsj” pointing of the 
payload instrument by allowing two rotations of the instrument relative to the 
carrier vehicle. For extensive Earth-pointing usage an alternative coarse 
gimbal arrangement [fig. 1(b)] may be preferable in which an elevation and a 
lateral gimbal provide simpler crosstrack capability when the space shuttle is 
in an inverted position relative to the local vertical. This coarse gimbal 
arrangement is assumed throughout this report. The second or vernier pointing 
assembly is made up of magnetic actuators for suspension and fine pointing, 
roll motor segments, and an instrument mounting plate around which a continuous 
annular rim (fig. 2) is attached which provides appropriate magnetic circuits for 
the actuators and the roll motor segments. The payload assembly, consisting of 
the payload instruments, the mounting surface, and the annular rim, is suspended 
in the elevation coarse gimbal without physical contact, and provides vernier 
attitude fine pointing and roll positioning of the instrument as well as six- 
degree-of-freedom isolation from carrier motion disturbances. In addition, 
the second assembly has a rim centering mode in which axial and radial rim 
position sensors located at each actuator station are used to center the rim 
axially and radially between actuator pole faces. This mode allows coarse gimbal 
slewing for retargeting, for Earth pointing, or for backup coarse gimbal 
pointing. 

Nominal operation of the ASP system for solar or stellar pointing first 
involves coarse gimbal pointing with coarse (wide field of view) sensors onboard 
the ASP system (or with .sensors onboard the carrier vehicle and relative gimbal 
angle information). The rim centering mode is activated during gimbal slewing. 
After coarse alinement, coarse roll positioning is accomplished by means of the 
rim roll motor and a relative roll sensor located on the rim together with a 
carrier sensor or a sensor onboard the ASP system. After coarse attitude aline- 
ment, verni jr fine pointing is initiated. For this mode, errors obtained from 
fine (narrow field of view) attitude sensors located on the ASP system, either 
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as part of the instrument or remote, are nulled by small magnetic-suspension- 
and fine-pointing-actuator torques applied to the annular rim. Also, 
translation rim centering is accomplished for this mode. 

Nominal operation of the ASP system for Earth pointing initially involves 
establishing the correct azimuth, elevation, and roll attitudes and attitude 
rates (slewing) in the rim centering mode. After appropriate smooth instrument 
slewing is established, vernier pointing (similar to that described for solar 
am. stellar pointing) is accomplished. 

The basic objectives of the present study are: 

• Derivation of a detail id mathematical model of the space shuttle/ ASP 
system. 

• Design of control laws in order to obtain optimum performance. 

• Prediction of the pointing performance when subjected to stochastic 
crew-motion and other disturbances, and sensor and actuator noise. 

In reference 1, a linear mathematical model was developed for the space 
shuttle/ASP system. This model ignored the payload- to-shuttle coupling. Sensor 
noise and actuator noise were also ignored. An azimuth-elevation coarse gimbal 
configuration was used, and the coarse gimbals were assumed to be fixed relative 
to the shuttle. However, as pointed out earlier, the "elevation- lateral" 
coarse gimbal arrangement is more suitable, especially for Earth-viewing 
e.xperiments , since it avoids gimbal lock. There'^ore the present study assumes 
the latter coarse gimbal arrangement. The coarse gimbals are no longer fixed, and 
control systems are designed for moving the coarse gimbals in such a way as to 
minimize the summed norm of the magnetic suspension centering errors at the four 
magnetic bearing stations. The roll freedom and the magnetic roll actuators* 
which were not modeled in reference 1, will be included in the present study. 
Sensor noise and actuator noise are also included. The mathematical model 
obtained in the present study describes the system in considerable details. All 
the internal couplings, cross products of inertia, etc., are included in the 
model. This semiannual report deals mostly with Task 1 described in the original 
research proposal, namely, linear analysis and design. 

Part 2 of this report contains a summary of research accomplishments and 
planned publications. Part 3 contains the derivation of a mathematical model 
using Lagrangian formulation. The design of controllers for the coarse gimbals 
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and the magnetic actuators is accomplished in part 4. Part 5 contains normalization 
of the equations of motion and inclusion of measurement noise and stochastic 
crew-motion model. Present and future research efforts are summarized in part 6. 
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1.1 Symbols 

A, Ap A2» a, Aj, A 2 


A2 


a 



>4, A. 

Bj, Bji 

B2 

B. 

J 

b 


Cl 

Cp (X) 


AB 


E 


E 

o 

e. 

J 


F. 

3 


F 


c 


F 


f 

c 


System coefficient matrixes 

Crew-motion filter coefficient matrix 

Points on the payload assembly annular rim which 
correspond to the magnetic bearing stations on 
lateral coarse gimbal 

Magnetic force coefficient matrixes 

Controller coefficient (rate gain) 

Coefficient matrixes 

System input matrixes 

Crew-motion 

Magnetic bearing stations on lateral coarse gimbal 
Controller coefficient (proportional gain) 

-S ‘'0,) 

Cross-product matrix of vector x. [x x y = Cp (x) yj 

Transformation matrix from A-coordinate system 
into B-coordinate system 

Angular velocity transformation matrix defined 
in equation (2) 

Normalized parameter matrix for the payload 

Column vector containing all zeros except a 
"1" in the jth position 
T 

(Fxj’ ^y j ’ forces at the jth 

bearing station 

Crew-motion force vector 

(with other subscripts) generalized forces 
corresponding to the subscript coordinates 

Normalized crew-motion force vector 
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gl* g2 

H , H , H 
Pi P 2 S oq 


I , I , I , etc. 

XX yy xy 


K , K 
gl g2 


L , L 
gl g2 


m 


m 


m , m 
gl g2 

m 

P 

“1 

0 


0 . 

1 


Elevation and lateral coarse gimbals 
Transformation matrixes 

(with subscripts) Inertia matrix defined in 
equation (4) 

(without subscript) Identity matrix 

Elements of Inertia matrix 

Objective function 

Coarse gimbal rate gains 

Coarse gimbal proportional gains 

m m + m 
s gl g2 

Mass of the shuttle 

Masses of coarse gimbals 

Mass of payload 

mm /(m + m ) 

P P 

Center of mass of payload 

Origin of the coordinate system (X^^, Y^, Z^) 
Center of mass of the shuttle 


*^1 

P2 

Q 

q 

‘Ir 

r 

T 

Tl, Tj, T3 


T , T , T , T 
s' gl g2 P 


gl 

z ) 

g2 


Origin of the coordinate system (X^ , , 

Origin of the coordinate system (X^ , Y^^, 
Center of instrument mounting surface 
Vector defined in equation (48) 

Generalized coordinate 

Radius of instrument mounting surface 

Transformation matrix 

Transformation matrices for 1, 2, 3 Eulerian 
rotations 

Kinetic energies of shuttle, elevation gimbal, 
lateral gimbal and payload respectively 
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con 


ABC 


ID 


(Xj, Y.. 2.) 

(X , Y , 2 ) 

^ s s s 

(X , Y , Z ) 
81 81 81 

(X , Y , Z ) 
82 82 82 

(X , Y , Z ) 

P P P 


(X , Y , Z ) 

q q q 


a , a 

p q 




e 

X 

V 

P 


Crew-motion torque vector 

Shuttle altitude control torque vector 

Normalized crew-motion torque vector 

White-noise source with unit power spectral density 

Noise covariance intensity matrix 

Linear velocity of point "A” relative to B-coordinate 
frame, expressed along C- coordinates 

White noise vector 

Measurement noise 

State vector 

Inertial coordinate system 

Inertial coordinate system 

Coordinate system fixed to shuttle 

Coordinate system fixed to gimbal gj 

Coordinate system fixed to gimbal g2 

Coordinate system fixed to payload, through its 
center of mass 

Coordinate system fixed to payload through Q 
Shuttle altitude vector 
Payload altitude vectors 
Normalized shuttle altitude vector 
Error variable 

(three-axis) centering errors at the jth 
magnetic bearing station 

Payload assembly position error 

Controller output 

Normalized payload assembly position error 
(with subscripts) damping ratios 
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'ss 


m 


ni2 


(♦ 3 . e^. 


(4> , B , \l> ) 

1- P P 

^ q q q 


♦ 1 

♦ 1 

e 

g 


g 


0' 


'g 


gl 


82 


’AB 


Subscripts 


i 


gl 

82 

P 

q 

s 

o 


Covariance matrix of x 
Steady-state value of ^ 

Coarse gimbal driving torques 
Components of 

Elevation coarse gimbal angle 
Components of 
Components of 
Target pointing angle 
Normalized value of ^ 

g 

Lateral coarse gimbal angle 
Target pointing angle 
Normalized value of 6 

g 

•Angular velocity of g frame relative to s-frame 

Angular velocity of g 2 frame relative to gj frame 
Inertial angular velocities, natural frequencies 
Positio.i of point A in coordinate system B 


Inertial frame (X., Y., Z.) 

V 

Inertial frame (X. , Y. , Z. ) 

11 11 11 

Elevation coarse gimbal 
Lateral coarse gimbal 
Payload 

Payload frame (X , Y , Z ) 
Shuttle 

Norminal or measured parameters 
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2. SUMMARY OF RESEARCH ACCOMPLISHMENTS 

2.1 Research Accomplishments 

The system configuration considered in this study is shown in figure 3. 

A detailed mathematical model of the space shuttle/co.-rse gimbals/fine pointing 
system was developed using the Lagrangian approach. The model automatically 
includes all the internal couplings. The equations were linearized around a 
set ot target angles in order to maintain the mathematical tractability of the 
problem, and to facilitate the derivation of appropriate control laws. Control 
systems were designed for the elevation and lateral coarse gimbals, which 
minimize the sum of the norms of the 5-dimensional centering errors at the four- 
bearing stations. Thus the coarse gimbals automatically and continuously follow 
the payload in the solar-, stellar-. Earth-pointing modes. Control laws were 
designed for the magnetic suspension and fine pointing actuators, which 
minimize the sum of the norms of the 5-dimensional magnetic forces at each 
bearing station. The canonical mathematical model for the space shuttle/ASP 
system has 11 degrees of freedom (order 22). Sensor noise and actuator noise 
arc also included in the model. This model, when coupled with the crew-motion 
model (ref. 2' of order 24, results in the total system order of 46. A computer 
program has been prepared to generate the system matrixes, and to solve the 
the resulting covariance equation. 

2.2 Planned Publications 

A paper on the mathematical model of the Annular Suspension and Pointing 
System and Predicted Performance Under StocK'.stic Crew Motion Disturbances 
and Sensor Noise is being planned, for submission probably to the AIAA 
Journal of Spacecraft and Rockets. This paper shall include forthcoming 
numerical results. 


3. MVTHEMATICAL MODEL DEVELOPNENT 


3.1 Mathematical Model Development 

(X., Y., Z.) is an inertial system centered at 0. . (X , Y , Z ) is a 
^111 1 s s s 

shuttle-fixed system, centered at 0^, the shuttle center of mass (fig. 3). 

0 nominally coincides with 0.. (X , Y , Z ) is a coordinate system 
s 1 gl gl gl 
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fixed to the elevation coarse gimbal, centered at Pi, and is obtained by 

one rotation 0 about X . The axis X remains parallel to X . The 
g gl gl 

coordinate system (X , Y , 2 ) is fixed to the lateral coarse gimbal. 

g2 g2 g2 

This system is centered at ? 2 , and is obtained from the gj system by one 

rotation 0 about Y . All the coordinate systems will hereafter be 
g g2 

referred to by their subscripts, e.g. i, s, gi , g 2 , etc. 


Another inertial coordinate system (X. , Y. , 2. ) is also centered at 
0^, and is obtained by two rotations ((n, 0i (X, Y sequence) from the 
i-system. and 0; represent the target pointing angles relative to the 

i-frame. 


The payload assembly consists of a cylindrical payload instrument, the 
instrument mounting surface, and the annular rim (fig. 4). (X^, Y^, 2^) is 

a coordinate system fixed to the payload assembly, centered at 0, its center 
of mass, ''h-.- payload-fixed coordinate system (X^, Y^, 2^) is parallel to 


(X , Y , 
p p r 
mounting surface. 


and is centered at point Q, which is the center of the instrument 


The coordinate svstem (X , Y , 2 ) is obtained fiom (X., Y. , 2.) by 

' s s s 1 1 1 


0 


in X, Y, 2 sequence. The vector 


a set of three rotations <Ji , 

•y S ^ ^ 

a = (c , e , C ) thus defines the shuttle attitude. T.ne pavload svstem 

(X , Y , 2 ) is obtained from the ii svstem bv a set of three rotations 
P P P 1 . . ^ 

(X, Y, 2 sequence) , 0 , . Thus the vector a = (4i , 0 , V ) defines 

‘ ^ p’ p p p p p p 

the payload attitude. 


In the material that follows, denotes the coordinates of point A 

in the B coordinate system; v.__ denotes the linear velocity of point A 
relative to the B-frame, expressed along the C-coordinates , The vector u> 
denotes the inertial anguiai velocity of the body expressed along the body- 
coordinates, while denotes relative angular velocity vector (as given in 

the list of symbols). The matrix D denotes the transformation matrix from 

A-coordinate system to B-coordinate system. 


5.2 Summary of Mathematical Model 


The mathematical model is obtained via Lagrangian formulation. The 
complete system consists of two separate systems: i) shuttle and coarse 

gimbals, and ii) payload assembly. Since these two systems are physically 
separated by magnetic actuator gaps, it was advantageous to derive separate 
models for each system, and to interconnect them via magnetic actuator forces. 
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The voarse pointing assembly consists of three objects; shuttle, elevation 
gimbal, and lateral gimbal. Expressions are derived for the kinetic energies 
of the three objects in terms of the inertial velocities of their centers of 
mass, and their angular velocities relative to the inertial frame, expressed 
along the body coordinates of each object. The eight degrees of freedom are: 


^0 ii ~ ^^0 i’ ^0 i* ^0 i^ ’ 
s s s s 


a , (fi and 
s’ g 


The gimbal bearings are assumed •'■o be frictionless, 
the six degrees of freedom are: 

a and = (X_. , Y„. , Z_. 

p Oil ^ Oil’ Oil Oil'' 

Linearization is performed about: 


0 . 

g 

For the payload assembly. 


^0 i ■ “s ^ *^g " '*’1’ 

^Oii " ^Oii 


0 

g 


® 1 » 



0 , 


Total degrees of freedom are 14; however, since only the relative position of the 
payload assembly, relative to the coarse pointing assembly is of interest, 
three degrees of freedom are redundant. Thus there are elev'en minimal degrees 
of freedom. 


3.3 Kinetic Energy of the Shuttle 

The inertial angular velocity of the shuttle resolved along the s-system is 

( 1 ) 


where 


0) = E a 
s s s 



'cos 0 

cos 

sin 


s 


s s 

E = 

-cos 0 

sin it cos il) 

s 

s 


s s 


_sin 0 


0 


s 



The linear velocity of 

0 

s 

in the i- frame 


0 

0 

1 J 


ii ~ ^0 i* 
s s 


The kinetic energy of the shuttle is given by; 

IT 1 ® T T * 

T = i m v^^ . . v» . . + = a E I E o 
s 2s OiiOii 2s s sss 
s s 


( 2 ) 


( 3 ) 
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where 


's = -I 


XX 

s 

"^XY 

s 

X 

l-H 

1 

XY 

s 

Iyy 

s 

-hz 

X2 

"^YZ 

'zz 


3.4 Kinetic Energy of the Elevation Coarse Gimbal 

The inertial angular velocity of the gimbal gi resolved along 
gl frame is 


where 


u) = n + D 0) 
g gl Sgl s 


= (*V 0 , 


D = T, (♦ ) 

sgl ^ g' 


= angular velocity of elevation coarse gimbal w. r. t. the shuttle) 
g 1 

The linear inertial velocity of the center of mass (located at Pj) expressed 
along i-coordinates is 

T 

V ..=v ..+D. (coxc 1 (6) 

^Piii Ogii is •- s ^pisJ ^ ^ 

= v^ . . + D. ^ H„ E a (7) 

0 11 IS Pi s s 
s '■ 


where 


H = -C (5 ) 

Pi P Pls" 

Kinetic energy 

.^1 T 1 T _ 

T ==m V ..v.. + -(D I 0 ) 

gl 2 gl PI 11 Pill 2 gl gl gl 

where I is defined similar to (4). 
gl 

3.5 Kinetic Energy of the Lateral Coarse Gimbal 

The inertial angular velocity of lateral gimbal resolved along the 
g 2 ~axes is 

0) = ■♦■1) U) 

g2 g2 glS2 gl 


u)=n+D Jl+D D 0) 
g2 82 8182 81 8182 sgl s 
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where 

o J 

fi = (0, 9 , 0) = angular velocity of lateral gimbal (along g£ 

coordinates) relative to g^ frame. 

The linear velocity of the center of mass (P 2 ) relative to s-frame 
resolved along gx is 

V =v +v 

P2Sgi Pisgi gi ^P2gi P2glgl 


V = D V 
P 2 SS Sgi P2Sgi 


V =V_. + V +<*>xr 

P2IS OglS P2SS S ^ 2 ^ 

therefore 

T T 

V.-V„.+D (v +n xr )+ujx(r +D ? ) 

P 2 IS O^lS Sgi PiSgi gi T2gl S ^PiS S gi P2gi 

Point Pi is fixed relative to the shuttle; 

V =0. 

PlSgl 

Since 


C = (0, 0, z ) , 
P2gl P2gl 


after simplification. 


where 


V =D. v-..+y+H E a 

P 2 IS IS 0 11 ' P 2 S s s 


Y = (0, -sin 4 , cos (b ) z 
^ ^g ^g P2gl 


H = -C (c + y) 

P 2 S p "'»pis 


Kinetic energy 


^ 1 T IT. 

T =-m V .v.+iU) I 0) 
g2 2 g2 P2 IS P 2 IS 2 g2 g2 g2 


3.6 Kinetic Energy of Payload Assembly 


For the payload assembly, it was found to be advantageous to use the 
coordinate system (X^, Y^, Z^) for obtaining the equations of motion. The 
(X^, Y^, Z^) system, which is payload fixed, is obtained from the ii -system 
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via three rotations (((» , 0 , ) = a in 1, 2, 3, sequence, and by translating 

M M M M 

the origin to Q which is the center of the instrument mounting surface. The 
vector then defines the payload attitude. The payload angular velocity 
relative to the inertial frame, resolved along the q-axes, is 


u = E a 

q q q 

where is defined in a manner similar to (equation 2). 


(15) 


The kinetic energy is given by 

^ I T 1 T ^ 

T = ^m_v_. v^. +su I u) + m 

p 2 p Q iq Qiq 2 q q q p 


|v„^. H 0) 1 (16) 

L Q iq oq q J 


where 


H = -C (C ) (17) 

oq p oq '■ ^ 

and = inertia matrix of the payload assembly about the q-axes. 

3 .7 Generalized Forces 

The generalized forces are obtained by increasing one generalized coordinate 
at a time, from q^ to q^ + 6^^, all other coordinates being fixed, and by 
examining the expressions for work done. Thus the generalized forces for the 
generalized coordinates are as follows: 


Coordinate (s) 

Generalized 

force 



T 

F - D D 

T 

4 


0 1 
s 

c sgi gl g2 

j=i^ 




T 


"0 



mi 



a 

T + T 

0 

- D ^ 

T 

s 

c con 


s gl 

m2 



.0 . 


.0 . 


4 

- r 

X 

T 

(D D 

F 



s gl gl 

g2 


g 


Qil 




z ^ 

j=l 

4 


o- (D F.) 
Bjgl gl g2 J 


T ~ ^2* y* Co • ^ ^ ■ 

n>2 ^ Bjg2 J 

4 

j = l ^ 


Z ^ F- 

jti J 
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(C^ = crew location in the s-system) 

3.8 Equations of Motion 

The total kinetic energy of the shuttle and the coarse pointing assembly is 

T = T + T + T 
s gi g2 

1 o<ro 1“TT ® 

2 s 0^ ^ 0^1 2 s s s s s 

10 •T' O'PO 'F O 

+ = m^ + D. H E a )* (r + D. H E o ) 

2 gi 0^^ is PI s S'' 0^^ 1 s PI s s'* 

+ 5 Cn + E a I (fi + D E o ) 

2 gi sgi s S'' gi '• gi sgi s s'' 


+ 5 ra (D. C . + Y + H E a (D. ; • + y + H E a ) 

2 g2 IS 0^1 P2S s s"^ ^ is 0^1 P2S s s'^ 


1 T 

+ = ^ ^ +D D w) I (fi +D n 

2 g2 glg2 gl glg2 Sgi S g2 g? gig2 gl 


+ D D 0) ) 
glg2 sgi S'* 

The equations of motion are given by: 


(18) 


d_(^ T\ _ ^ 

dt \ ^J " 3 q^ " *"q^ 


(19) 


where q , r=l, 2, n, are the generalized coordinates and F are the 

generalized forces. Application of (19) to the kinetic energy in equation (18) 
results in the equations of motion, wnich are nonlinear and highly complex. 
Linearization is performed about the equilibrium point.. 

^ = 0, Og = 0, (p^ = (Pi , 

and the following equations are obtained: 

mCn +(m H +m H )a+hz m (p = F . (20) 

0^1 gl Pi g2 P2S" s P2gl g2 ^g 

where 


m = m m ♦ m 

s gl g2 
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I 

I 


t 


(I+m H^H +171 


+ D 


T . 

D__ ♦ D,, I D.. ) a 


S gl Pi Pi g2 P2 S P2S S gi gi Sgi 111 g2 111 

+ (m H ^ + m H ^ . + D ^ I ei ♦ m z H ^ h 

gl Pi g2 P2 s' 0^1 » gl gj 1 g2 p2gi P2 S 


T •• T •• 

+ D. . I D ei) $ + D. . I e? 9 = F 

1 il g2 glg2 g 1 11 g2 g a. 


( 21 ) 


(m z h^) C . + (ei^ I D +2 H m 

g2 P2gl Og P2gl P2S g2 

T T •• T T 

+ 6] D I 0.. ) a + (Ivy + ei D I D ei 

^ glg2 g2 111 S XXgi glg2 g2 gig2 ^ 


+ m z‘- ) (f + Co I D ei 0 = F, 


g2 P2gl g 


g2 glg2 


g 4>g 


9 + I D ei + 02 *^ ID D a = F0 

YYg 2 g g2 glg2 g g 2 glg2 sgi s g 


( 22 ) 

(23) 


In equations (20) - (23), the symbols , a , , 6 denote their 

*Jgi s g g 

incremental values about the equilibrium point. D , D denote the 

^ sgi’ gig2 

corresponding transformation matrices evaluated at constant angles <l>i and 

01, and D. . = D D (evaluated at (|)i, 0i). 

^ 111 gig2 sgi 

For the payload assembly, the linearized equations around 
“q ^ “• 'Qii ^ « 

are : 


m (C^. + H *a ) « F^«. (24) 

p ^^Qii oq q' cQii 

M T •• 

I a + m H = F (25) 

q q p oq ^Qil aq 

•• 

substituting for from (24) into (25), 

(I-mH^H)a=F -H^ F^^ (26) 

q p oq oq' q aq oq tQii 

But the moment of inertia of the payload assembly about the p-axes is 

I = I - m H ^ H . 

P q p oq oq 

Also, for the linearized system. 


a 3 a 

q P 

Therefore, equation (26) becomes: 

I a = F - H ^ 
p p aq oq cQii 


(27) 
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Considering the two position equations, (20) for and (24) for 

, only the difference between the payload assembly position and the 
coarse pointing assembly is of interest. 


Define 


' Sail ■ Sii 


Before linearization, 

C . = D. . C . + D. . D J D ; + D.. dJ ; 

P2I1 111 0^1 111 IS s gi P2gi 111 1 s PIS 

After linearization of the second derivatives of this equation about the 
equilibrium point, simplification yields 

• ft T* ft# ft# ap •• 

. = D. . (e ♦ ) - H a - hz ♦ (29) 

^0^1 1 ii Qil'' P2S s p2gi "^g 

This value of Cq ^ is substituted into equations (20) - (23) to obtain the 
following: 


•• bX 

e ♦ D.. \ H 

111 \ m Pi 


m + a 
s gi 


H /a 
P2S / s 


/m, + m \ D. . 

( s gi j •• lil 

•D. . \ /hz - H a = — — F 

111 \ m / P2gi g oq q m c 


1 * — 1 E f- 

"■ V 1^1 " 


where is the crew-motion force. The remaining equations of motion are: 

(I+m H^H -m H^H +D D ♦ dJ. I D. . ) ? 

^ s gi Pi Pi gi Pi P2S sgi gi sgi 1 ii g2 111 s 

T T T - T 

+ (m H ‘ ♦ m H ‘ ) D. . e + (D^ „ ei 

gl Pi g2 P2 s 1 ii s gi gi 

T T •• T •• 

♦ D. . I D ei - m H hz ) (^ •♦■ D. . I ®2 

1 11 g2 gl82 gl Pi P2gl g 1 il 82 8 

- (m H ♦ m H /) D.^. H ^ = T^ + T 

81 Pi 82 P2S 1 ii oq p c con 


T cos 

m2 


T Sin ^1 
m2 ^ 


1 ■ S ^^Bjs^(^i ii *"j) 


17 


I 


(m H ^ + m H ) 
81 PI 82 P2 s 


82 P28l 
♦ (I, 


P 

). . c •*• 
1 11 

T 

+ ei D 


d/. y F. ♦ 

^ pi J 


(ei 

I 

D + ei 

• D 


81 

S8l 


T 



•• 

1 

I 

D ei) 

4 > 

81 82 

82 

8182 

8 

•* 


. T 


ei 0 

- m 

z h 

D. 

8 

82 P281 

1 


; X F 
c c 


I D.. ) a 

81 82 82 111 S 


q* 99 

♦62 I D eie-m z 
82 8182 8 82 P 

T ^ 

= T - ei T C ) D 

mi jti P 


Bj8l 81 82 j 


82 7 

z h 

mp P 281 


D. . y F. 


YYg 2 8 


e + 62 I D 


82 8182 
4 


•• T •• 

<j) + 62 I D. . a 

^ rt ^ rt ^ < 


82 111 s 


= t - 92 * Z C ) F ( 33 ) 

m2 P Bjg2 J 

I a = y C (;..-;) F. ( 34 ) 

P P f^l P oq^ 3 

Equations ( 30 ) - ( 34 ) d6fine th6 compl6t6 linearized motion of the system. 
Position Sensor Outputs : 

An axial and a radial position sensor is located at each magnetic actuator 

station. For computing the magnetic actuator forces, it is necessary to 

compute the centering errors at each actuator station. Since the payload 

assembly rim is continuous, it is necessary to compute the points on the rim 

directly under the actuator stations. Tliis can be done by minimizing over rim 

roll angle ijj, the distance between a point on the rim (assumed to be a circle) 

and a magnetic actuator station. For the present linear analysis, however, 

this can be achieved by a simpler method. The linearized radial centering 

errors are simply the radial components of c. The axial centering errors 

are obtained by finding the angle between the bearing plane and the rim plane. 

The vector • is transformed ("projected”) into the q'- coordinate 

system (with e = 0), q'- system being same as the q-system, but without the 

ij/p = rotation. The Z-component of the transformed vector plus the 

axial (Z-axis) components of translation error c, gives the total axial 

centering errors. Let 5 . (x) denote the X -axis centering error at the 

J 82 

jth actuator station, etc. 
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First ignoring the translation. 




T2(0p) Ti(^p) 


D. . 
Ill 




glS 2 


''Bjg2 


(35) 


After linearization and simplification, the following expressions for the 
sensor outputs are obtained: 


5 1 (x) = 6 3 (x) = (36) 

5 2 (y) = 5 4 (y) = Ey (37) 

(z) = e - r (cos (})i 9 + sin 4ii ) - r (9 - 6 ) (38) 

62 (z) = e + r (cos 9^ $ ♦ sin 9^ sin ijii 9 - sin cos il/ ) 

Z S 5 5 


+ r (Og cos 9i - (J)p) (39) 

63 (z) and 54 (z) are obtained from 5 ^ (z) and 6 2 (z) by changing the 
sign of r. 


4. CONTROL SYSTEMS DESIGN 


4.1 Coarse Gimbal Control Systems 

The coarse gimbal assembly, which consists of the elevation and lateral 

coarse gimbals, must automatically follow the payload assembly. This is 

accomplished by the use of two brushless DC torquers. The objective of the 

control systems for the torquers is to move the coarse gimbals in such a manner 

that the payload assembly is kept properly centered in the magnetic bearings. 

Thus an appropriate design approach is to obtain the command values of 

6 and 9 which will minimize 

E g 

4 

j = s (x) *63^ (x) * S 2 ^ (y) * 5 4 ^ (y) * 5 (x) (<o) 

j=i J 

This is done by substituting for 5^ from equations (36) to (39), and for 
c from linearized version of equation (28), and making 



3 J 
3 e 


g 


* 0 
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The solution gives the following command signals: 


♦ - 41 = i fz e - 

^gC ,2 ^ J.2 cos2 e, '■ ^281 y 

P2gl ^ 

- r2 cos 01 (cos 01 ♦ sin 0i sin 4 >i 0^ 

®gc *^ 5 ^ 

It can also be proved that 


cos 01 (cos 01 ♦ 


g 



- sin 01 cos 4 >i ( 41 a) 

( 41 b) 


1 


^ z2 + p2 cos2 01 ^2gl 

P2gl 


tp,., « 2 (y) - 


r cos 01 


|52 (O - 64 (z)[] 


(42) 


and 


0 - 0 = 
gc g 


5 1 (z) - 5 3 (z) 
2r 


(43) 


Equations ( 41 a), ( 41 b) are used for solution of the state equations, and 
( 42 ), f 45 ) are used in actual implementation. 


The torques to be generated by the coarse gimbal torquers are then given by 


O 



- K ^ 

♦ L 


- ♦ ) 

( 44 ) 

mi 

gl g 

g} 

gc 

g 


0 

♦ L 


<Z> 

1 

( 45 ) 

m2 

g2 g 

g2 

gc 

8 

where K , K 

gl g2 

are rate 

gains 

and 

Lgi, are proportional gains. 

The 


desired response can be (approximately) obtained by making 


K 

gl 


T T 

lyv D 

Xygl glg2 


I D ei 
g2 gl82 


2p u 
gl 81 



^XXgi 


T ^T 
ei D 


8182 


I D ei 
82 8182 



'vVg2 



u , and 
82 


L 

82 


^Y^'g2 



(46) 


( 47 ) 
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4.2 Magnetic Actuator Control Systems 


Magnetic forces acting in and directions (nominal) at each of 

the four bearing stations are used for controlling the motion of the payload 
assembly. Tims there are twelve control forces. Ignoring the rotation- 
translation coupling, the desired responses are obtained by making 


and 


where 


j=i 


XJ 


y F • 

^ yj 

j = l 
j = l 


c y* F. 

j^i ^ 


= ">10 


rf = -I X 

po ap 


f=(F -F,F -F.F 


+ F 


Z 4 


22 21 


F. ) 


23 X2 X 4 y 3 Xl 


mi = mm^Cm + m J , = -C^ (;^^). X^ = (X^^, X^^, X^^) 


cx' ey' ez 


X =2p u) e + w^e, etc. 

EX XXX XX 


The subscript "o" implies the measured or estimated value of the parameter. 

Thus the control force equation can be written in the following matrix-vector form: 


where 


Acf F “ q 


Letting 


F=(F ,F ,F ,F ,F ,F ,F ,F ,F ,F ,F ,F ) 

'■ Xl* X2* X3’ X4’ yi* y2* Y 3 Yk zi’ 22' 23* Z4-' 


q = r mio X*^ , - (I X C mio X )*^ 1 
M LIU e r po ap Co ^ ■* 

A^P = appropriate 6 * 12 coefficient matrix. 


c r 1 T , 

F = F = i mio ei X 
Xl X3 2 ^ e 

F = F = ^ mio 62 ^ X 
y2 y4 2 e 


(48) 


( 49 ) 
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The remaining equations are: 




or 


1 1 
0 0 


0 0 0 0 0 0 
1 1 0 0 0 0 
0 0 111 


0 0 0 -1 0 


1 -1 


1 0 -1 

0 0 0 


q - 


X 2 


X4 


yi 


X3 


22 


23 

, F 
I- Z14 


The solution which minimizes the norm of F is given by 
- ' T i ' "Ti”!" 

F “ (^CF^ f ^^CF^ I ^ 

which, after simplification, gives: 

F “ ^ 2 * ®3^ (I ^ "*10 

X 2 4r ^ po ap Co e 


F = -F 

X 4 X2 


F = F 

yi X 4 


F » F 
y3 X 2 

Fjj = 5 niio e,’’ ej’’ dpo \p * So S’ 

T 1 

m • A A A 1 _ 

e 2 r 


Z 2 4 J e 2 r ^ ^ po ap co c 


1 T IT 

F - z min 03 X - -:r— 62 (I X + C min X ) 
zi 4 ‘lu j j. 2 r ^ po ap co c' 


1 T IT 

F - z min 63 X + -:r- ei (I X ♦ C mjo X ) 
Z 4 4 iu J E 2r * ^ po ap Co e 
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This can be oppressed as: 
F » >4q 


F ♦ F 

XI X3 

F ♦ F 
Y2 Yk 

0 

0 

0 

0 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 
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( 


where 


.T 1 


1 

ei 0 

ei 0 0 e2 0 e2 

1 

^3 

1 

2^3 

1 

2^3 

2 

1 

[O -^C2 

1 1 1 
0 -263-2^3 0 ^ej 0 


-ei 

-e2 

ei 


~ * i+4 

eT 

i+8 


i = 1, 2, 3, 4 


F. = -4. i = 1, 2, 3, 4 (59] 

These are the expressions for desired magnetic actuator forces. 

4.3 Space Shuttle Attitude Control System 

The attitude of the shuttle is assumed to be controlled by such devices 
as CMG's or AMCD. The control torque is given by 

T *-(I+m H^H +m H^H ) \ (60] 

con " s gi Pi Pi g2 P2 s P 2 S'' a. '• ^ 


where 


s s s s 


\ = “♦ *S * “% *5- 

s s s s 


5. COMPUTATIONAL ASPECTS 


5.1 Equations of Motion UTien Actual Payload Parameters Differ From the 
Nominal Ones 

The control systems designed in the proceeding section are based on nominal 

or measured values of the parameters i; , n and I , which are, in general, 

oq P P 

different from the actual values. It is desirable to investigate the pointing 
performance when such errors exist. To keep the number of variables small, 
it is assumed that the same per unit error 6 exists between the actual and 
measured (or nominal) values of these parameters, i.e. 

«) nip . (1 ♦ 6) nip^, Ip - d * «) Ip„. 


23 



I 


Then equation ( 34 ) becomes 


•• 


or 


I a * (C - C ) min X - I X 
p p ^ ; ;o“^ G po ap 


(1 ♦ 6) I a = 6C^ mio X - I X 
^ po p CO e po ap 


which reduces to 


• • 

a = 
p 


1 


■; T I C, min X , -f- 

1 6 po Co ^ “ e 1 ♦ 6 ap 


( 63 ) 


( 64 ) 


( 65 ) 


Similarly, equation ( 30 ) becomes 


111 


e Ai2 o ♦ A13 A15 = — 


1+6 


(1 ♦ 


6m 

m ♦ m ' c 
po 


( 66 ) 


where A^2» Aj3 and A^ are easily identified by examining equation ( 30 ). 

In an attempt to normalize the equations of motion to some extent, the 
Euclidian norm I’QqQil* hereafter denoted by JcH, was selected as the 
normalizing parameter. The following variables and parameters are defined: 

Eq = miO Ipo'' UP. Cl = V = c/ld. Sj " 

♦ig = '’c ' Vi'l- 'c ' Vli'l 

The normalized equations of motion then become 

•e •• •• » 

V + Aip + Aj3 ♦ Ai 5 Op 


1 


= B 11 f - , " T (1 * 5 E^) X 

1+6^ m+mv 

po 

•••••• •• •• 

A21 V ♦ A22 ♦ ^23 4 iig ♦ A24 6lg ♦ A25 Op 


U 1 


T » t 
con c 


1 

ku 


mi 

T COS 

m2 

T sin ^1 
m2 


( 67 ) 


- y c (C.. ) D.^. F. 
p Bjs'' 1 1 ] 


T T T 

(m H ♦ m H * ) D. . . 

gl Pi g2 P2 1 -1 ^ 


(1 ♦ O nip^ II cH 


F . ♦ c * f 

^ j c c 

j-1 


( 68 ) 


•• •• •• •• •• 

A3I * A32 65 ♦ A33 (JUg ♦ A34 0 ig ♦ A35 Op 
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m 


'id L'”! 

ei 

4 

j=i 

•1 

z h’'D..T 

E 

^i] 

P2gl 

111 

j = l 

e« 

•0 



+ A43 <}»ig + 

A44 

61 

^g 


Z (c„.„ ) d'' 


F. - 


g2 


po 


(69) 


1 




u\2 


The matrices are given in the appendix. 

•® 6 1 

“ = T^~r E C, X - -7-^ ^ 

p 1+6 o^v 1+6 op 


T C ) F. (70) 

Cl j = i 


(71) 


Equations (67) - (70) contain the terms F^ /Id. in their right-hand 

sides. From equations (59) and (48) 




where 


Li = 


^3x3 


LCi 


^3X3 


E 


c| 


(72) 


(73) 


o r -* 

Thus the Fj/|U| terms in the right-hand rides of equations (67) - (70) can 
be normalized. The only quantities that remain to be normalized are x 


and T , which are given by equations (44), (45), (41a) and (41b). 
in^ 


mi 


After simplification, division by |c| yields 
<P - 4> 

1££ li = 

Icll 

' * P2gl 


r 

+ r^ cos^ 01 L 


2 V - r^ cos 01 ^(cos 01 (^1 


P2gi y 


!(■ 


g 


Ic| 


(cos 01 , sin 01 sin <J»i, - sin 0i cos (()i) (74) 


0 


6 - 6 0 

^ (0, cos (fii, sin 4»i) 6 

|c| Ul 


(75) 


Hence it is seen that iji and 0 , and therefore x , and x cannot be 

gc gc ml m2 

normalized in this fashion, and the parameter r or |c| must be specified 
in addition to the normalized parameters 
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Cl = C^q/|c|, r/|c| 


and 


(76) 


E„ = "lO Id' Ip„-‘ 


5.2 Representation of Measurement Noise 

Control systems for the coarse gimbals process the measurements of 

e o 

6 . , j = 1, 2, 3, 4, , 6 , and generate appropriate torque signals. Control 

3 S 8 

systems for magnetic actuators process the measurements of 6., j = 1, 2, 3, 4, 
o o 3 

(for E and e) and a^, in order to generate appropriate voltages for the 

magnetic actuator coils (force equations 50 - 57). Space shuttle attitude 

control system makes use of the shuttle attitude and attitude rate measurements 

0 

(a and a ). All these measurements contain certain amounts of noise. The 
payload assembly attitude and rate are estimated by a star-tracker 

system fixed to the payload which possibly makes use of a Kalman filter. The 
star-tracker model currently available is discrete-time^. However, as a first 

O 

attempt at the problem, continuous measurements of a^, are assumed to be 
available, and are represented by adding continuous white noise to the actual 
quantities. The covariance of the additive noise is to be made equal to that 
reported in reference 3 using the "optimum system." The results obtained 
using this approximation are expected to be conservative (showing worse perfor- 
mance than actual) because the actual estimation error is not white, but is band- 
limited. All other measurements are contaminated with continuous white noise, 
and appear in all control terms, e.g., X^^, etc. Sensor dynamics have 

not been included in this first trial, but appropriate representations will be 
included later. A nonlinear model for the coarse gimbal torquers has been 
supplied by Sperry Flight Systems, and will be used in task 2, the latter part 
of the grant. The linear part of this model contains the representation 
of torque as being proportional to the motor current. Since the control loop 
contains current feedback, the lag due to inductance can be theoretically made 
as small as desired. For the second-order model currently available, the 
eigenvalues represent dynamics which are much faster than the rest of the 
system; therefore the torque outputs may be assumed to be instantaneous in 
this first attempt at the solution. 

Additi'w measurement noise contaminates the sensor outputs at the 
four suspension gaps (6j, j = 1, 2, 3, 4), coarse gimbal angle rates, 
payload attitude and attitude rates. The coarse gimbal angular rates are 
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measured by tachometers, and the relative displacement rate (e) is generated 

by using a lead circuit. Thus the measurement noise, represented as white 

noise processes, appears in all the control terms; e.g., X , X , etc. Each 

V p s 

noise term is normalized by division by ||c||. 

The crew-motion forces and moments for each crew activity can be described 
by passing a unity power spectral density white noise through a filter of order 
four (or less) for each force and each moment.^ The complete equations of 
motion can be symbolically written as 

A x\ = Ai Xi + A2 Xi + Bi f^^ + Di w^ (77) 


f = T Xo 
cr ^ 


(78) 


X2 = A2 X2 B2 u 


where 


Xi 





T.T 

“p^ 


f 

cr 



t 

c ^ 


(79) 


X2 = state vector of crew-motion filters 


u = white noise with power spectral density = l/||cp 


Let 


w^ = measurement noise vector (normalized) 


Xi = Yi 


Xi = Y2 = Yi 


O 



= Y2 


Y 2 = A"1 Y 2 + A-1 ^2 Yi 


+ A-1 


Bj T X 2 + A'^ 


Di 


w 

m 


X2 — A2 X2 82 u 


Thus, defining 


T o T T X 

X = (Yi‘, Yi‘, X2')‘ 


(80) 

(81) 

(82) 
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equations (80), (81) and (82) can be written in the form 

X = A X + I V (83) 

where v is a vector white noise process. 

The resulting covariance equation for this system is 

O 

H = AJ3 + J^A^+BVB^ (84) 

where 

E [v v^] = V 6 (t - t) (85) 

This covariance equation can then be solved for ^ 55 > steady 

state value of ^ , by the method described in reference 4. 

6.0 PRESENT AND FUTURE EFFORTS 

A detailed linearized mathematical model has been derived for the space 
shuttle/ASP system via Lagrangian formulation. The model incorporates all 
internal couplings, and sensor noises. The elevation-lateral coarse gimbal 
configuration has been assumed. Control laws have been designed for the coarse 
gimbals, which minimize the summed norm of the suspension gap errors at the 
bearing stations. Thus the coarse gimbals have the capability of automatically 
and continuously following the payload in order to point at a target which is 
moving with respect to the space shuttle. Minimum norm control laws have 
been designed for magnetic suspension and fine-pointing actuators. A computer 
program has been written for generating the system matrices, and for solving the 
resulting covariance equation. 

In the immediate future, preliminary computational results will be 

obtained for this model using the program mentioned above. Computation will 

be performed using data for two or three t>'pical payloads, and statistical 

pointing accuracies will be obtained as functions of errors in estimating the 

payload parameters, and sensor noise variances. The star-tracker models which 

are currently available will be reviewed, and a suitable model will be chosen 

for inserting in the present mathematical model. (The present mathematical 

model assumes that the payload attitude estimation error is continuous white 

noise, with covariance reported in reference 3. This is a conservative 

assumption and should give worse- than-actual pointing accuracies.) A detailed 

star-tracker/payload attitude estimator model would be discrete-time; therefore 

its interconnection with the continuous shuttlc/ASPS model will have to be 

investigated. One solution would be to discretize the continuous system, and 
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to use equivalent discrete white measurement noises and crew-motion filter input 
noise. A simpler and computationally less expensive solution would be to 
approximate the discrete-time star-tracker/estimator system by a continuous 
system. 

The above computational results would also yield the effect of the payload 
motion on the shuttle. If this is found to be significant, a model for a second 
payload/ASPS/coarse gimbal system will be added to the present model, and the 
computation will be repeated. The two payloads are connected only through the 
shuttle motion; therefore these efforts need be undertaken only if the payload 
motion is found to significantly affect the shuttle. 

If the knowledge of flexible modes of the shuttle is available, this shall 
be included in the present model as the next step. The linear analysis shall 
be complete with this step. 

The next task is to perform nonlinear analysis. The nonlinearities to be 
considered include: 1) coarse gimbal torquer nonlinearities due to ripple, 

cogging and friction; 2) magnetic suspension and fine-pointing actuator 
nonlinearities: square-law, deadcone, hysteresis, coil switching, etc.; 

3) sensor deadcones; and 4) shuttle attitude control system vernier jets, which 
operate in on-off fashion. The nonlinear system will be represented by a linear 
part and time-invariant nonlinear blocks. Computational results will be 
obtained via Monte-Carlo simulations. Finally, time permitting, an analysis 
of earth-pointing performance shall be performed using the linear model. One 
way to perform this could be to evaluate the ASP's local pointing performance 
when tracking a relatively fast-moving target. 
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APPENDIX 


Coefficient Matrices: 


( m m + m \ 

-SJ- h - -5 H_ ) 

m Pj m P2S/ 

/ID + m \ 

Ai3 = - D. . I— SJ-) hz 

iM\ / P2gl 


"'^15 “ ■ (1 5 ) 

D. . 

B,- = ^ 


Ml 


'21 
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= (m H ^ * m H )d^. . 

\ gl Pi g2 P2 s / 111 


+ m H - I D <■ D*^. . I D. . 

s gl Pi Pi gl Pi P 2 S sgi gl sgi 111 g2 111 




A23 = ( I + D^. . I 

' sgi gl 111 g2 


D \ 

glS2 / 


ei - m Hn hz_ 

gl Pi P2gl 


= D . . I eo 
111 g2 ^ 

Aoc = - (m H 




Ag 1 = m z_ _ . . 

g2 P2gl 111 


h" d". . 


''^32 - ®1 


(l D + I D.. ) 

\ gl sgi gig2 g2 111/ 
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A33 = ^YVn ■*■ ei D ID Cl 

^Xgi ^ glg2 g2 gig2 * 

T 

A34 = eo I D ei 
^ ^ g2 glgZ ^ 

A35 = - m 2_ D^. . Cl (1 + 6) 

g2 P2gl 111 
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A42 ~ 6? I D. . 

^ g2 111 

T 

A|. 3 = 60 I D ©1 

g2 glg2 * 


A44 = I 


YYg2 
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Figure 1. Annular Suspension and Pointing (ASP) System 
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Figure 2. Vernier Pointing Assembly Configuration 
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Figure 4. Payload Assembly 
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